Models of organometallic complexes for optoelectronic applications 



O 
(N 



o 

o 

• I— I 
Oh. 



> 

On 
00 
(N 

P 

o 
o 



X 



A. C. Jacko, Ross H. McKenzie, and B. J. Poweljj 

Centre for Organic Photonics and Electronics, School of Mathematics and Physics, The University of Queensland 

Organometallic complexes have potential applications as the optically active components of or- 
ganic light emitting diodes (OLEDs) and organic photovoltaics (OPV). Development of more ef- 
fective complexes may be aided by understanding their excited state properties. Here we discuss 
two key theoretical approaches to investigate these complexes: first principles atomistic models and 
effective Hamiltonian models. We review applications of these methods, such as, determining the 
nature of the emitting state, predicting the fraction of injected charges that form triplet excitations, 
and explaining the sensitivity of device performance to small changes in the molecular structure of 
the organometallic complexes. 

PACS numbers: 



INTRODUCTION 

The development of organic light emitting diodes [1| 
and dye-sensitized solar cells [2] two decades ago has lead 
to an explosion of interest in the optoelectronic proper- 
ties of organic materials for use in photovoltaic and elec- 
troluminescent applications. [sUsl For example, a material 
with the right set of optoelectronic properties could revo- 
lutionize solar energy production, replacing fragile, high- 
temperature processed silicon with durable, solution pro- 
cessed, ink-jet printable plastic electronics. ^6n9|] However, 
the efficiency of optoelectronic devices needs to increase 
in order for them to become commercially viable beyond 
niche applications. Improving their efficiency while re- 
taining the desired chemical and physical properties has 
proven challenging. [9|, llO| 

Some of the most promising candidate materials for 
light-emitting devices are organometallic complexes. Fig. 
[T] shows the structures and spectra of two typical 
organometallic complexes used in OLEDs. The presence 
of a heavy transition metal ion core allows triplet states 
to radiatively decay (via a large spin-orbit coupling) . It is 
widely supposed that upon injection triplets and singlet 
are formed in the ratio 3:1. Hence phosphorescence al- 
lows for a more efficient utilization of injected charges, [llj 

The technologically useful properties of organometallic 
complexes, namely their absorption, emission and charge 
transfer properties, are determined by the excited states 
of the complex. For example, the process of electrolumi- 
nescence involves the injection of positive and negative 
charges into a bulk sample of the complex, which re- 
sults in the ionized complexes. These charges migrate 
until they combine on a single complex, leaving it in an 
electronically excited state. This excited complex then 
decays back to it's ground state by emitting a photon. 
A key design criterion is that the radiative decay must 
be much faster than any non-radiative decay processes. 
In OPV devices, incoming photons create excited states, 
i.e., an electron- hole pair. These charges must then be 
separated to perform useful work. Thus, in OPV devices. 



one wants charge-seperation to be the preferred decay 
pathway for excited states. Fig. [2] shows a schematic of 
charge recombination in an OLED. To understand and 
optimize this process (and the reverse process in OPVs), 
one must first understand the various excited states in- 
volved. Fig. [3] shows a schematic of the likely decay 
pathway in typical phosphorescent OLED materials. 

A key property of the excited states of organometallic 
complexes is their character, i.e. the degree of singlet or 
triplet character and the degree of ligand centered (LC) 
or metal-to-ligand charge transfer (MLCT) character. By 
understanding which parts of the complex influence each 
excited state, one might begin to engineer the spectrum 
of a complex by tuning the chemistry of the complex. 
Notice that there are two bands in both of the absorp- 
tion spectra shown in Fig. [TJ The higher energy, more 
absorptive band is identified with LC (tt — tt*) transitions. 
This higher energy (usually UV) band is present in spec- 
tra of isolated ligands, and thus its identification as a LC 
transition is uncontentious. The lower energy, lower ab- 
sorptivity band is usually associated with MLCT singlets 
as it is seen only in complexes. p/&] The classification of 
the emitting state has attracted considerable interest and 
debate, and it has been labeled variously as MLCT, LC 
or MLCT-LC hybrid. [Ts] The emitting state is believed 
to be a triplet, but the precise degree of triplet character 
has long been a question of interest. [19] 

The brightness of the MLCT singlet is due to the hy- 
bridization of (occupied) metal orbitals and (unoccupied) 
ligand tt* (LUMO) orbital. This mixing delocalizes the 
excitation over the whole complex, allowing the excited 
state to have a transition to the ground state of sim- 
ilar brightness to a localized tt* ^ tt (LC) transition. 
This shift of electron density from the metal to the lig- 
and is known as tt back-bonding. This interaction has 
several experimental signatures, including an increase in 
the length of the C-C bond in the ligand and a softening 
of the ligand C-C stretch frequency as well as the bright 
MLCT transition, [io^ii 

The presence of a heavy transition metal core in an 
organometallic complex produces a strong spin-orbit cou- 
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FIG. 1: Absorption and emission spectra of Pd(2-thpy)2(top 
panel) from Ref. [l2, and ppyPt(dpm) (bottom panel) from 
Ref. 13. In both panels the dashed line shows the absorp- 
tion spectra and the solid line shows emission spectra. The 
brighter higher energy absorption bands are usually identified 
as LC (tt -^ TT*) singlets, the less intense lower energy band 
is usually identified as MLCT singlets. The emission is typi- 
cally considered to be predominantly triplet, and, depending 
on the complex may be LC, MLCT or a hybrid of the two. 
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FIG. 2: A schematic of the process associated with the oper- 
ation of a LED. From top to bottom: The neutral complexes 
(C) are in their ground states. An electron is injected from 
the right, and a hole from the left. These ionic complexes are 
now in their new (charged) ground states. The electron and 
hole combine on a single complex (by either electron transfer 
between the LUMOs of adjacent complexes, or hole transfer 
between the HOMO-metal hybrid orbitals), forming a singlet 
or a triplet excited state (C*), while the other site returns to 
the neutral ground state. Finally, the excited state relaxes to 
the neutral ground state, either by emitting a photon or by 
non-radiative mechanisms. 
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pling, since the magnitude of the spin-orbit interaction 
increases proportional to Z^ (where Z is the atomic num- 
ber of the element). [17| The large spin-orbit interaction 
in the presence of heavy metal ion causes a hybridization 
of singlets and triplets. [19] This has several important ef- 
fects: First, the states dominated by singlets and triplets 
are now coupled by a rapid intersystem crossing (~50 fs) 
via their MLCT components. [15, 17] Second, spin-orbit 
coupling allows triplet-dominated states to decay radia- 
tively. This allows for phosphorescent OLEDs, increasing 
the excitation-to-photon conversion efficiency by harvest- 
ing triplet excitations. [11] 

Finkenzeller et a/.[23[ have argued that the degree of 
zero-field splitting in the emitting triplet reflects the 
amount of MLCT character in the state. They further 
state that there is 'an empirical correlation between the 
amount of ZFS and the compound's potential for its use 
as emitter material in an OLED\ They conclude that 
increasing the MLCT character of the emitting triplet is 
'a necessary but not sufficient condition' for improving a 



FIG. 3: A schematic of the decay pathways in a phosphores- 
cent material based on the results of Ref. [ij- Injected charges 
form either an MLCT singlet excited state i^MLCT) or a hy- 
brid LC-MLCT triplet state ^LC+^MLCT). MLCT singlets 
undergo a fast (fs) intersystem crossing (ISC) to the triplet 
state, the speed of which is determined by the magnitude of 
the spin-orbit coup ling and the degree of MLCT character 
in the triplet state, p^-flq The triplet excited state relaxes 
to the ground state (GS), either radiatively, with a rate kn, 
or non-radiatively, with a rate kNR- The photoluminescent 
quantum yield (PLQY) of the complex is then ^ 



kR + kNf^ • 



complexes suitability for use in OLED devices. [23| 

Understanding the character of the excited states may 
help make better materials for organic electronic devices. 
Furthermore, it may allow us to understand why empir- 
ical correlations, such as those of Finkenzeller et a/., are 
true and when we can expect them to be valid. Further, 
there remain many open questions about the properties 
of organometallic complexes including: 



• Why are their emission properties so sensitive to 
small chemical changes, cf., e.g.^ Ref. 'S"? 

• What is the character of the emitting triplet state 
- to what extent is it MLCT, or LC? 

• What is the non-radiative decay path of the emit- 
ting state? What part do unoccupied non-bonding 
metal orbitals play in non-radiative decay? 

• What is the mechanism for ultrafast singlet-triplet 
MLCT transitions? Is there a conical intersection 
associated with this intersystem crossing? 

• Are the excited states localized over all the ligands, 
or localized on one? How does this change with 
perturbations to the symmetry of the complex? 

• What fraction of injected charges result in triplet 
excitations? 

Clearly, theoretical input has an important role to play 
in answering these, and other, questions. To date, ma- 
jority of modelling of organometallic complexes has been 
atomistic first principles quantum chemistry. However, 
another approach is to model Hamiltonians that capture 
the main features of broad classes of materials. In the 
remainder of this paper we compare and contrast how 
these two complimentary approaches have been used to 
address some of the above questions. 



FIRST PRINCIPLES QUANTUM CHEMISTRY 

First principles approaches have the advantage of be- 
ing materials specific. This can be invaluable for under- 
standing the role of chemical changes in the complexes. 
However, it is also a disadvantage in that it requires 
many different complexes to be studied before one can 
hope to access these trends, much like experiment. The 
organometallic complexes used in OPV and OLED appli- 
cations present a number of challenges to atomistic mod- 
elling. They are quite large, which limits the use of high 
quality ab initio methods - although some calculations 
have been reported. [2J Yet electronic correlations do 
play an important role, particularly in determining the 
excitation energies that determine the optical properties 
of the complex. Ideally one would like to be able to pre- 
dict excitation energies to within chemical accuracy {i.e., 
within ~ ksT). By far the most widespread approaches 
to the first principles modelling of organometallic com- 
plexes have been density functional theory (DFT) for 
ground state properties and time-dependent density func- 
tional theory (TDDFT) for excited states properties (see, 
for example, Refs. I5[ l25l432l ). (TD)DFT gives reasonable 
accuracy at a low computational cost. While these meth- 
ods do not provide chemical accuracy for the excited state 
properties they have allowed for trends across related 



complexes to be studied. In particular, many researches 
have attempted to identify the character of excitations, 
i.e., whether they are ligand centered (LC), metal-to- 
ligand charge transfer (MLCT) or inter-ligand excita- 
tions, with (TD)DFT (see, for example, Refs. SlUI-ii). 

While TDDFT is a useful resource, current implemen- 
tations have some well known shortcomings. |33l438J For 
example, it is well known that current (TD)DFT imple- 
mentations do not correctly capture electronic correla- 
tions or describe charge transfer processes well.|37H39| 
Current functionals can successfully predict semiconduc- 
tor band gaps of inorganic crystals to within 0.1 eV, but 
much of this accuracy can be attributed to a fortuitous 
'cancellation of errors'. [3J In the context of organometal- 
lic complexes, one of the worrying consequences of these 
problems is that (TD)DFT tends to overestimate the de- 
localisation of electron orbitals. |37l l39| In order to bet- 
ter understand both the attributes and the limitations 
of (TD)DFT calculations we finish this section by briefly 
reviewing three calculations [5|, l25|, l26[ - highlighting on 
both the insights gained and the limitations of the calcu- 
lations. 

Hay [25] used TDDFT to investigate the effects of per- 
turbing an Ir(ppy)3 complex by changing one ligand, 
replacing one phenylpyridine (ppy) with either acetoy- 
lacetonate (acac) or benzyolacetonate (bza). He found 
that the HOMOs of each of these complexes have around 
50% metal character, the rest being ligand tt, while the 
LUMOs were predicted to be purely ligand tt* orbitals. 
Thus, he argued that excitations in these complexes are 
hybrid MLCT-LC. When one of the ppy ligands is re- 
placed by acac or bza Hay predicted only small changes in 
the HOMO energy. However the HOMO-1 and HOMO-2, 
which are degenerate in Ir(ppy)3 due to its C3 symmetry, 
are split and stabilized by the symmetry perturbation 
introduced by the new ligand. In Ir(ppy)3 the LUMO 
was predicted to be a combination of all three ligands, 
as were the LUMO+1 and LUMO+2: a degenerate Eg 
pair. In Ir(ppy)2(acac) the LUMO and LUMO+1 were 
predicted to be degenerate and localised on the two re- 
maining ppy ligands. This degeneracy is rather surpris- 
ing as this molecule has C2 symmetry and thus no sym- 
metry required degenerices. In contrast the LUMO+2 
was predicted to be locahsed on the acac ligand. Hay 
found that in Ir (ppy) 2 (bza) the LUMO is localised on 
the bza ligand. He argued that this localisation of the 
LUMO in Ir (ppy) 2 (bza) is the reason that this com- 
plex has a much lower PLQY than Ir (ppy) 2 (acac) or 
Ir(ppy)3 (< 0.01 in Ir(ppy)2(bza),[40] compared to 0.3 
and > 0.4 in Ir(ppy)2(acac)[40] and Ir(ppy)3[41], respec- 
tively). However, it is not known whether this change 
comes via a supression of the radiative decay rate or an 
enhancement of the non-radiative decay rate. 

Obara et al. [26] investigated a series of iridium com- 
plexes with both bi- and tri-dentate ligands. They used 
TDDFT calculations to investigate the large transition 



dipole moment from states that are predominantly triplet 
found in many organometallic complexes. Following pre- 
vious discussions of rhodium complexes, |42i |43| they con- 
sidered two mechanisms by which spin-orbit coupling 
might allow triplet radiative emission. One is direct cou- 
pling between a LC triplet and MLCT singlet. The other 
is the coupling between the MLCT component of a pre- 
dominantly LC triplet and the MLCT singlet. They then 
found that the direct LC to MLCT spin-orbit coupling is 
unlikely to be large enough to explain the bright triplets 
observed experimentally. Obara et at. conclude that 
a MLCT triplet to MLCT singlet coupling is the most 
reasonable explanation for the intense phosphorescence 
of the complexes they investigated. These results sup- 
port the hypothesis of Finkenzeller et at. that increasing 
the MLCT character of the triplet will increase its ra- 
diative decay rate and hence its potential as an emitter 
material. [26l 

Lo et al. (5[ investigated a series of iridium complexes 
with aryltriazolyl-based ligands experimentally and with 
DFT calculations. They found, experimentally, that 
small chemical changes cause large variations in both the 
radiative decay rate (discussed further in the following 
section) and the non-radiative decay rate. They sug- 
gested a correlation between the measured photolumines- 
cent quantum yield (PLQY) and the calculated molecular 
orbital energy gaps.[5| However, they found that this cor- 
relation only remains as long as the non-radiative decay 
rate does not vary much. In addition, they found that 
the HOMO of Ir(ppy)3 has five times as much electron 
density on the phenyl ring as on the pyridyl ring, while 
the LUMO has three times as much electron density on 
the pyridyl ring as on the phenyl. Thus, substituents 
that stabilize the phenyl ring will lower the energy of the 
HOMO more than that of the LUMO. Similarly, desta- 
bilizing the pyridyl ring raises the energy of the LUMO 
more than that of the HOMO. Therefore both stabilizing 
the phenyl ring and destabilizing the pyridyl ring result 
in a blue-shift of the spectrum. Thus, they were able 
to explain why one can tune the HOMO and LUMO of 
Ir(ppy)3 somewhat independently. 



MODEL EFFECTIVE HAMILTONIANS 

Another approach to modeling the optoelectronic 
properties of organometallic complexes is to construct 
a model with fewer states but an accurate treatment 
of the electronic correlations. This contrasts with first 
principles calculations, which include several basis states 
for each atom but neglects some electronic correlations. 
The small number of degrees of freedom in such semi- 
empirical models allows one to make fewer approxima- 
tions on the interactions and correlations in the model. 
It also allows one to identify key trends that describe 
broad classes of materials. This approach has proven 



itself incredibly powerful in wide areas of materials sci- 
ence. For example, the Anderson single impurity model 
can describe a wide range of systems including magnetic 
impurities in metals, quantum dots in semiconductor 
heterostructures, carbon nanotubes, and single molecule 
transistors. (41145! 

In principal an effective model Hamiltonian is found 
by starting with the exact Hamiltonian and 'integrating 
out' high energy states. [46| This procedure is computa- 
tionally expensive, |47| so often one simply chooses a re- 
duced basis set, motivated by the physical processes one 
wishes to capture. [46| DFT can be used to estimate the 
values (or trends in values) of some of the parameters of 
these effective models. The model Hamiltonian can then 
solved, retaining correlations that the approximate DFT 
functional does not include. |48M5 2] 

Identifying the frontier orbitals which dominate the 
photophysics is one of the most significant steps of the ef- 
fective model approach. In this reduced basis set one can 
define an effective Hamiltonian with just a few parame- 
ters. Conjugated polyenes have been investigated in this 
way via the Hiickel, Hubbard, Heisenberg and Pariser- 
Parr-Pople models. [46] This approach has been applied 
to organometallic complexes, for example mixed valence 
binuclear systems including magnetic atoms in proteins 
(Hubbard and double exchange models) (Ref. |53|), molec- 
ular magnets (Ref. 54), Anderson impurity models for 
cobalt based valence tautomers (Ref. 1551 ). and a series 
of metal-cored bipyridine complexes (Ref. 'lOl). It has 
also been shown recently that this approach naturally ex- 
plains the sensitivity of the photophysical properties of 
organometallic complexes to small chemical changes. [14 1 

To correctly describe the character of the excited states 
the model must capture the key interactions. There are 
many important features of the system that might be 
included in such a model, for example electronic 'hop- 
ping' terms between the frontier molecular orbitals, di- 
rect Coloumb interactions between electrons in those or- 
bitals, spin interactions, and relativistic effects such as 
spin-orbit coupling. The relative energy scales of these 
various interactions will define the composition of the ex- 
cited states and therefore their properties. Again, some 
examples will help to illustrate both the power and the 
limitations of this approach. 

Kober and Meyer [19] studied a simple model Hamilto- 
nian for P3 symmetric bipyridine (bpy) complexes. Their 
basis consists of 3 HOMOs and 3 LUMOs (two of E sym- 
metry and one of A). The HOMOs are of mixed metal- 
ligand character, while the LUMOs are pure ligand tt*. 
The parameters of the model are the HOMO-LUMO gap, 
the splitting of the A and E HOMOs, the splitting of 
the A and E LUMOs, the singlet-triplet splitting and 
the degree of spin-orbit coupling. To make progress this 
model has to be fitted experimental spectra. This is a 
significant limitation. However, the fit reproduces the 
experimental data quite well. On the basis of this model. 
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FIG. 4: Experimental phosphorescent radiative rates versus 
calculated singlet-triplet gaps scaled by atomic number Z, as 
in Ref. ITq , whose data we show here. Ref. [ly discussed 
this correlation for the iridium and platinum complexes, all 
of which have similar values of Z. One can easily extend this 
to complexes with very different values of Z (and therefore 
spin-orbit coupling) by including the Z^ dependence in the 
transition dipole moment. [17,] This shows that reductions in 
the singlet-triplet gap are correlated with increases in the ra- 
diative decay rate. The solid red line is a fit of the data to the 
predicted kn oc Z^ / /S^E'^j, behaviour. The dashed blue line is 
a line of best fit to kn oc {Z^/AEItV with a = 2/3, indicat- 
ing that although that, over this data set, the predictions are 
only followed qualitatively. 



they concluded that although there is appreciable singlet- 
triplet mixing, labelling states as 'singlet' or 'triplet' is 
still reasonable, and can be helpful in understanding the 
photophysics. They also found that although the initial 
excitation occurs to a symmetric combination of bpy lig- 
ands, the relaxed excited state does not necessarily have 
this symmetry. 

Haneder et al. [16] investigated the correlation between 
the Si-Ti gap, AEst^ and the radiative decay rate of 
the triplet, kR. They considered a three level system 
consisting of the So, Si and Ti states and calculated the 
radiative decay rate of the triplet state by including the 
spin orbit coupling perturbatively. They found (within 
their small sample of molecules) that there was indeed 
a qualitative trend similar to their prediction that kR oc 
AEgj.^ cf. Fig. m The strength of the spin-orbit coupling 
varies with the atomic number, Z, of the transition metal 
as Z^, while the radiative decay rate is proportional to 
the square of the transition dipole moment. Therefore, 
one should also expect that kR oc Z^, Fig. |4] shows that 
this is followed as to the same accuracy Haneder et aVs 
original prediction is. Thus, this very simple model leads 
to two straightforward design principles: (i) minimize the 
singlet-triplet gap to maximize the triplet radiative rate 
and (ii) maximise the atomic number of the transition 
metal. 

We have recently applied a semi-empirical approach to 



modelling the optoelectronic properties of organometal- 
Hc complexes. [14] We employed a basis of effective frag- 
ment orbit als, one pair of frontier ligand orbit als and one 
metal orbital. The Hamiltonian for this model explicitly 
included the hopping between the ligand orbitals and the 
metal, an effective spin exchange interaction between the 
ligands HOMO and LUMO, and the renormalised direct 
Coloumb interactions on and between every site. It is 
important to realise that these effective renormalised pa- 
rameters will in general be very different from the bare 
parameters used in first principles methods ^, ,48,] . Dif- 
ferences between complexes are encapsulated into the ef- 
fective parameters of this three orbital model. Neverthe- 
less, this model Hamiltonian can be solved exactly. 

For the parameters characteristic of the organometallic 
complexes used in OPV and OLED devices, the lowest 
excited singlet eigenstates are pure MLCT and LC, while 
the triplet states are typically strongly mixed. The pa- 
rameter that controls the character of the triplet states 
is the ^MLCT-^LC energy gap (as shown in Fig. [5]) 



J/4 



(1) 



where £* is the renormalized effective ligand HOMO- 
metal energy gap and J is the exchange interaction be- 
tween an electron in the ligand HOMO and an electron 
in the ligand LUMO. J determines the gap between the 
lowest singlet triplet excited states in the isolated com- 
plex. The ^MLCT and ^LC states are not eigenstates. 
Thus the size and sign of the gap between them deter- 
mines the character of the lowest excited triplet state, Ti. 
Note that a molecular orbital theory would predict that 
singlets and triplets are degenerate (J = 0), a prediction 
that is not borne out experimentally. 

Most of the parameters of our model can be extracted 
from either basic theory or from comparison to experi- 
ments on the isolated ligand [1^. However, the param- 
eter £* cannot be predicted from studies of fragments of 
the complex, it is a property of the complex, the partic- 
ular combination of metal and ligand. In particular it 
depends on the complex-specific properties the HOMO- 
metal energy gap e and the strength of the LUMO-metal 
Coloumb interactions. Fig. [5] shows the two important 
parameter regimes of the model. In one (s*/J < 1/4) the 
lowest excited triplet is predominantly LC, while in the 
other (e*/ J > 1/4) the lowest excited triplet is MLCT. In 
the latter case, the energy gap between the lowest excited 
triplet and lowest excited singlet becomes small. 

This variation in triplet state character can have pro- 
found effects on the optoelectronic properties of the com- 
plex. £* / J can have a dramatic effect on the triplet's 
radiative decay rate, as illustrated in Fig. [6l[14] As e* / J 
increases, the lowest triplet state's MLCT character in- 
creases and it gets closer in energy to the MLCT singlet 
state, both of which serve to increase the effects of spin- 
orbit coupling. In the regime in which the spin orbit 
coupling can be included perturbatively {e* / J ^ 1/4), 
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FIG. 5: Correlation diagram, showing triplet basis states on 
either side and the resulting hybrid orbit als in the center. The 
upper figure corresponds to the situation s* /J < 1/4, where 
the Ti state is predominantly LC, while the lower figure shows 
s* /J > 1/4 where the Ti state is predominantly MLCT. In 
the intermediate case e* /J :^ 1/4, both the Ti and T2 triplet 
states are strongly mixed MLCT-LC hybrids. 



the triplet radiative rate increases by more than an or- 
der of magnitude when e* /J is varied from 0.1 to 0.25. 
In this regime the radiative rate can be shown to vary 
with the singlet-triplet gap to the fourth power, which, 
on the scale in the figure, appears linear on a semilog- 
arithmic plot. [14] In the regime £*/J < 1/4, the degree 
of MLCT character in the triplet is perturbative, and we 
find agreement with the hypothesis of Finkenzeller et al: 
increasing the MLCT character increases the radiative 
rate of the emitting triplet. [23| 

We also determined the effect of the parameter £:*/J 
on triplet production probability following charge injec- 
tion. The characteristic time for singlet-triplet intersys- 
tem crossing is very short compared to the singlet and 
triplet lifetimes, so there will be a Boltzmann distribu- 
tion of excited states. The lowest singlet and triplet are 
well separated (compared to ksT) from other excited 
states, so one only needs to consider these two lowest 
excited states. When the singlet- triplet gap is much less 
than ksT we have the expected 75% triplet probability to 
25% singlet probability. On the other hand, if the singlet- 
triplet gap is much larger than /c^T, we have nearly 100% 
triplet probability. This means that injected charges will 
almost always decay to the ground state via the lowest 
triplet state if the singlet-triplet gap is larger than ksT. 
In our model a large singlet-triplet gap corresponds to 
the regime £*/J < 1/4, as shown in Fig. [71 
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FIG. 6: Triplet radiative rate k^ (lower curves, blue) and sin- 
glet radiative rate (upper curve, red) as a function of £*/J 
for various values of Hso solved to first order in perturba- 
tion theory. We have chosen {{r)iMLCT^ — {r)o) — 20 A to 
reproduce a singlet radiative rate of ^ 10* s~^. The triplet 
radiative rate increases exponentially as £*/J increases, up 
until the point where the lowest singlet and triplet are nearly 
degenerate £*/J ^ 0.25 (i.e. both MLCT) at which point 
the perturbative solution becomes invalid. As the strength 
of the spin-orbit coupling increases, the triplet radiative rate 
rapidly increases. Around £*/J ^ 1/4 perturbation theory 
breaks down, and a more complete treatment is required. 



Therefore we have seen that this model gives a firm 
mathematical basis to some of the key ideas in the lit- 
erature. The hypothesis of Finkenzeller et al. that the 
triplet is predominantly LC with a perturbative MLCT 
contribution is reproduced in the £*/J < 0.25 regime of 
the model. In this regime, increasing £*/J simultane- 
ously decreases the singlet-triplet gap and increases the 
triplets MLCT character. Both of these lead to an in- 
creased triplet radiative rate, in agreement with the hy- 
potheses in Refs. 23 and 16. Note that increasing £:*/J 
does not guarantee an increase in PLQY, as we did not 
attempt to predict changes in the non-radiative decay 
rate (as this requires a knowledge of the non-radiative 
decay path on the excited state potential energy surface). 
Furthermore, one should expect large changes in the non- 
radiative rate to be caused by small changes in the chem- 
istry as non-radiative decay is an activated process and 
so will depend exponentially on energy gaps. 



CONCLUSIONS 

First principles electronic structure and effective model 
Hamiltonians provide complementary methods to un- 
derstand the optical properties of organometallic com- 



— J=2.12eV 

— J=1.54eV 

— J=1.17eV 
— ■ J=0.88eV 




FIG. 7: Triplet probabilities versus the ^LC-^MLCT energy 
gap £* — J/4 for various values of J, corresponding to the 
ligands ppy (J = 2.12 eV), thpy (J = 1.54 eV), fluorene 
(J = 1.17 eV) and bzq (J = 0.88 eV), all at 300K. 



plexes with applications in OPV and OLED technolo- 
gies. (TD)DFT calculations have for a long time been 
the workhorse in the field and have led to important in- 
sights, such as, those discussed above concerning how 
changes to the symmetry or chemistry of a complex can 
effect its PLQY. However, semi-empirical models can also 
provide significant insights: particularly for understand- 
ing trends across broad classes of materials. Here we have 
highlighted just a few examples: The radiative decay rate 
is sensitive to small changes in the ligand chemistry be- 
cause small chemical changes drive large changes in the 
degree of MLCT character in the triplet state. The same 
physics determines the fraction of injected charges that 
form triplet excitations. 
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